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In this Letter, weak turbulence theory is used to investigate interactions among Alfven waves and fast and 
slow magnetosonic waves in collisionless low-p plasmas. The wave kinetic equations are derived from the 
equations of magnetohydrodynamics, and extra terms are then added to model collisionless damping. These 
equations are used to provide a quantitative description of a variety of nonlinear processes, including "parallel" 
and "perpendicular" energy cascade, energy transfer between wave types, "phase mixing," and the generation 
of back-scattered Alfven waves. 
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PACS numbers: 52.35.Bj,52.35.Ra,95.30.Qd,96.60.Pb,96.60.Rd 

Turbulence at length scales smaller than the collisional 
mean free path Xmfp plays a central role in a wide range of 
astrophysical and laboratory plasmas. In general, the analysis 
of waves and turbulence at scales < A-mfp requires the use of 
kinetic theory. However, in some cases fluid models are ap- 
proximately valid even at such collisionless scales. For exam- 
ple, if p = Snp/B^ <C 1, where p is the pressure and B is the 
magnetic field, then magnetohydrodynamics (MHD) provides 
an approximately correct description of the fast magnetosonic 
wave ("fast wave") when X < A-mfp and CO <C £i,, where A, is 
the wavelength and i2, is the proton cyclotron frequency. [[l|] 
Similarly, MHD accurately describes both Alfven waves and 
anisotropic Alfven-wave turbulence when n <C A, < A-mfp and 
CO ^ Hi, where n is the proton gyroradius. HI |3] MHD is 
approximately accurate in these cases because the dynamics 
are governed primarily by magnetic forces and inertia, while 
the pressure tensor and collisionless damping play only a mi- 
nor role. In this Letter, MHD is used to model turbulence at 
length scales ^ n and < Amfp and frequencies ^ i2, in low-p 
plasmas. To account for the strong collisionless damping of 
slow magnetosonic waves and the weak collisionless damp- 
ing of fast waves, |1] extra damping terms are added to the 
equations for the wave power spectra. Although this approach 
is only an approximation to the full kinetic behavior of the 
plasma, the comparative simplicity of MHD makes it possible 
to describe the physics within the MHD model in great detail 
and thereby gain useful insight into the full problem. 

The basic phenomenology of MHD turbulence depends on 
whether the turbulence is weak or strong, which in turn de- 
pends on the value of COi^T^, where cOjt is the linear wave fre- 
quency at wave vector k and Zk is the time scale on which the 
fluctuations at wave vector k evolve due to nonlinearities. If 
N/th*: ^ 1. then the turbulence is weak, the fluctuations can 
be approximated as a collection of small-amplitude waves, 
and the interactions between waves can be analyzed using per- 
turbation theory. |l3t|4t] On the other hand, if |co<:|t<: < 1, then 
the fluctuations are not wave-like and the turbulence is strong. 
In MHD, 111 is at least as large as ~ (A;5vj-)^', where Svj. is the 
rms amplitude of the velocity fluctuation at scale ' . Thus, 
the condition |cOjt|T<^ ^ 1 is satisfied provided |coa.| 3> k^vi^. 

An important point is that the weak and strong turbulence 
limits can apply to different components of the turbulence 
within a single plasma. S H 01 For Alfven waves, cOa = 



±k^VA, where va = Bo/\/4%po is the Alfven speed, po is the 
background density, and Bo = Bqz is the background magnetic 
field. As a result, Alfven-wave turbulence is strong for suffi- 
ciently small jfcj^l/fcx^ where k± = k k,z. On the other hand, 
Alfven waves with \k,\ > k± and dv^ <C va are weakly turbu- 
lent. Similarly, fast waves satisfy (x)k — ±kvA in low-p plas- 
mas, and are thus weakly turbulent provided bv^ ^ va- This 
Letter focuses on weak turbulence, but a method to account 
for strong-Alfven-wave turbulence is also described. 
The equations of ideal MHD are 
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where p is the density and v is the velocity. The specific en- 
tropy [o<: ln{pp^^)] is taken to be a constant (where y is the 
ratio of specific heats). Each fluid quantity is taken to be the 
sum of a uniform background value plus a small-amplitude 
fluctuation: B = Bqz + 8B, p = po + bp, p = po + 8p, and 
V = Vo + 5v, with v'o <C 8v ^ \'a- The (spatial) Fourier trans- 
forms of 5v and b = bB/ ^/47tpo can be written as 

Vk = Va,kea.k + Vf,kki_ + v^±z, 



bk = ba,kea,k + bf^kef,k, 

where ea.k = z x k_i, k±_ = k_i/kj_, and e/^k = ea.k x k/k. The 
Alfven-wave amplitudes at wave vector k are 

+ I , 

The fast and slow-wave amplitudes, fj^ and i^, are given by 
w =M u, 

where w = {fk Jk^4^ h)^ « = ihk,Vf,k,v,,k,bf,k), hk = 
CsPk/po, Cs = vTPo/Po is the sound speed, and pk is the 
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Fourier transform of 5p. The matrix M is an infinite series 
in powers of e = Cs/va — ■\/yP/2. To order e^. 
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where /j = cos 9, ri = sin 9, and 9 is the angle between k and z. 
The Fourier transforms of equations ([T]) through (O, expressed 



in terms of sf,af, and /r^, become 
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where the right-hand sides are the nonlinear terms, co^^ 



ikyVA, and, to lowest order in e. 
±kvA. 

The power spectra are 
(4(4)*) = Sfb{k-kr), 



defined by the equations 
±(4)*) = A±5(* - k,). 



and (/+(/+)*) - 
ensemble average, 
the energy per unit 
waves) propagating 
proportional to the 
waves propagating 
about the z axis 



Fi;8{k — ki), where (...) denotes an 
The quantity (S^) is proportional to 
volume in A:-space of Alfven waves (slow 
I in the ±z direction. The quantity fj. is 
energy per unit volume in fc-space of fast 
in the k direction. Cylindrical symmetry 
is assumed, so that 5^ — S^{k^,k^,t), 
and — F{kj^,k^,t). 
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In the weak-turbulence limit, the wave kinetic equations 
can be obtained from equations (01 through (|6]l using the stan- 
dard techniques of I3i, These equations express dS^ /dt, 
dA^ /dt, and dFk/dt as series in powers of e. As written be- 
low, the lowest-order terms in these series are °^ e^^, con- 
tain S^, and are associated with the slow-wave density fluc- 
tuation. Pi- ~ fe,.Vj. |(.po/co^^,, which is a factor e^'csc9 larger 
than the fast-wave density fluctuation, pjt ~ ^±i'/,*:Po/K)^i., 
when v^^k = Although proportional to e^^, these terms 

may nevertheless be small, because strong collisionless damp- 
ing 1 1] makes much smaller than and F^. In this Letter, 
the S-^'^Sf terms are retained, but the nonlinear terms contain- 
ing 5* at higher order in e are dropped, with the exception of 
and co^^ = the d{qy) term in equation O, which is retained for reasons 
discussed below. Of the terms that do not contain 5*, only the 
leading-order terms e") are kept. The wave kinetic equa- 
tions then become 
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^Jd^pd\d{k-p-q)l^d{q,)^{k^n7n)\ (A+-A+) + d{k, + p, + q)k,A^_pk {kA^F^,, + p,F^^A+ + qApA+ 
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+ (S+p + S-p) \h{q - k,)m^ {F, - A+) + 5(^ + k,)m^ - A+) + 5(^,0^2 (A+ - A+) 

+ dik, + q,W {A, -At) ] + §(<?, + k^k^At^^ (^,A, ) I - 2y+^+ 

= ^Jd^pd\d{k-p-q)S^9smH[?>{k^p-q)kqFp{F,^Fk) + d{k + p ~ q)k{kF^pF^+ pF,Fi,~ qF^pF^) 
+ d{k - p, + q,)kAkpc, (M+A,7 - p-A^jFk + qzA+Fk) + b{k - p, - q)kMkpc, {kA+F^ - p.FyFk - qA+Fk) 
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whemAkpc = A{k,p^q)=k'^{kJ+2p±^m+2q±^n)^,Mkpq = 
M{k.,p,q) — k^^[k_J + px(cosa— l)m + ksinan]^, a is the 
angle between z and q, F^g — F{qj_,~q^,t), and y^j^, j^j^, and 
Y/jt are the hnear damping rates. The partial derivative dF^/dq 
is taken at constant a, and the partial derivative dA^/dq, is 
taken at constant ^x- In the triangle with sides of lengths k^^, 
p±, and q±, the interior angles opposite the sides of length k±, 
p±, and q± are denoted Ok, Op, and a^, and / = cosOjt, m ~ 
cosGp, n = cosa^, 1 = sinOA^, m = sino^, and n = sino^. The 
equation for dA^^ /dt is obtained by setting Aj 
F^k, 5* ^ S^, and y^^. ^ Y^.A: in equation 

The "collision integrals" on the right-hand sides of 
equations Q through (|9]l represent the effects of resonant 
three-wave interactions and sum over all wavenumber 
triads involving k that satisfy the resonance conditions 
k = p + q and co^ = (Op + (Oq, where cOa is the frequency at 
wavenumber k. When = at some wave vector ki, the 
only non-vanishing terms in dA^ /dt are non-negative at ki. 
Analogous statements hold for A^, and F/^. Equations O 
through (|9]l thus ensure that the spectra remain non-negative. 
When the linear damping terms are dropped, equations d?) 
through (|9|l conserve the energy per unit mass / d^k (A^ + 



2Fk + St 



S 



J/2 



+2coseFi + e 



and the pseudo-momentum 
1(5+-5,)]/(2va). When 



Jd'k [Ai 

the equation for d{v^/dt is taken into account, it can be 
shown that resonant three-wave interactions also conserve the 
cross helicity (v-B) and momentum (pv,). 

At k, = 0, the only nonzero term in the collision integral in 
equation ^ is the term proportional to ?>{qz)- This term rep- 
resents interactions between three Alfven waves ("AAA in- 
teractions"), which transfer Alfven-wave ener gy at all k^ to 
lai-ger k±^ but not towards lai-ger jgj [loi ITTl iH [ll H In 
AAA interactions, each Alfven wave type (fl+ or a^) is cas- 
caded by the other Alfven wave type. Thus, if A^{k±,Q) — 
[where A^{k±,Q) denotes Aj evaluated at k^ = 0], then the 
AAA term in dAf /dt vanishes. A Zakharov transformation 

can be used to show that A'^{k j_,0) <^ k^" is a steady-state 
solution to equation (|8]l for k^ = in the absence of dissipa- 
tion, provided «+ = 6, as in the incompressible case. lll2ll 
When dissipation is included, these power laws become ap- 
proximate solutions for A* {k± , 0) within the inertial range. If 
A+ {k± , 0) ~ A ^ (k \ , 0) at the (perpendicular) dissipation scale 
("pinning" IH [isl]) and A+{kj_,0) > A' {kj_, 0) in the inertial 
range, then > 3 > for the inertial-range spectra. The 
Alfven-wave spectra at k- = are not affected by the value of 
A^ at nonzero k^ or by the slow-wave or fast-wave spectra. 

At k- = 0, the only nonzero term on the right-hand side of 
equation (Q is the term §(q'z), which represents the mix- 
ing of slow waves by Alfven waves, which transfers slow- 
wave energy to larger k± but not to larger This term 
is identical to the expression describing the mixing of a pas- 
sive scalar by weak Alfven-wave turbulence, with Sf replac- 
ing the passive-scalar spectrum. In the "imbalanced" case 
in which A+(A;x,0) ^ A^ {k±,Q) within the inertial range, 
the quantity (A+ +A^) in this "passive-scalar mixing term" 



can be approximated as simply A+. If A^{k±,0) °^ k^" , 
a Zakharov transformation can then be used to show that 
steady-state solution to equation (|7]l 
at k- = in the absence of dissipation. Thus, the slow-wave 
spectrum at k- = (and hence also the spectrum of a pas- 
sive scalar) mimics the spectrum of the minority Alfven-wave 
type, A^(A;x,0). Although all other terms in the wave kinetic 
equations containing at orders higher than e^^ have been 
discarded, the d{q,)Sf term in equation d?) has been retained 
because it can dominate as k- 0, since the other nonlin- 
ear terms and the linear (Landau) damping term vanish in this 
limit. [Because strong coUisionless damping keeps Sf small 
at nonzero k^, the cascade of slow-wave energy to larger \kz\ 
arising from interactions among slow waves is neglected in 
equation dTji.] 

The term °^ e^^5(pz) in equation (O represents "phase- 
mixing." Slow-wave density fluctuations at k- = cause the 
Alfven speed to vary in the directions perpendicular to Bq- As 
a result, Alfven-wave phase fronts travel at different speeds 
on different field lines, transferring Alfven-wave energy to 
larger k±. [16] (Density fluctuations at k^ — associated 
with passive-scalar entropy waves would have the same ef- 
fect.) Phase mixing and AAA interactions both cause a 
perpendicular cascade of Alfven-wave energy. The relative 
strength of these two processes varies with 9. For example, 
if e"25±(/tx,0) ~ A+(A:x,0) ~ A"(A;x,0), then phase mixing 
dominates the perpendicular cascade when \kz\ ^ k± while 
AAA interactions dominate when ^x ^ l^zl- 

In equation (|9]l, the terms proportional to 9 sin^ 9 represent 
interactions between three fast waves ("FFF interactions"). 
The FFF terms are the same as the collision integral for weak 
acoustic turbulence up to an overall multiplicative factor 
proportional to sin^ 9. As sin 9 — > 0, the acoustic-like FFF 
interactions weaken because the fast waves become less com- 
pressive. [17] Energy is transferred from small k to large k by 
FFF interactions. [6, 17] The resonance conditions for FFF in- 
teractions require that p and q be parallel or anti-parallel to k, 
indicating that FFF interactions transfer energy along radial 
lines in A:-space. lli, 17| The terms containing A^p^ in equa- 
tions (|8]l and (|9]l represent interactions between two Alfven 
waves and one fast wave ("AAF interactions"). The terms 
containing M^pg in equations (O and (|9]l represent interactions 
between one Alfven wave and two fast waves ("AFF interac- 
tions"). Whenfex ^ \kz\, AFF interactions cause A *(A;x, A;,) to 
become approximately equal to F{k±,±\kz\). [17]. The com- 
bination of FFF and AFF interactions results in a "parallel cas- 
cade," i.e., a transfer of Alfven-wave and fast-wave energy to 
larger |^,|. iH 

The E^^5p terms in equation (|9]l represent the "resonant 
scattering" of fast waves into either new fast waves or Alfven 
waves of equal frequency but different wavenumber. The 
e^^b{k — q) term in equation (|9]l acts to isotropize F^. The 
e^^Sp terms in equation dH) other than the "phase mixing" 
term are also denoted "resonant scattering" terms, and repre- 
sent the conversion of an Alfven wave into a new Alfven wave 
or fast wave of equal frequency. If 5^ > A^ and Sf > Fi^, 
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then resonant scattering and phase mixing are the most rapid 
nonhnear processes in the p ^ Hmit. |8] On the other hand, 
in colHsionless systems, Landau damping can reduce S"^ suf- 
ficiently that resonant-scattering is weak. (Phase mixing in- 
volves at A;, — where Landau damping vanishes and thus 
can be very efficient even in colHsionless systems.) 

The "resonant-scattering" term in equation (O °^ £^^5(fc- + 
qz){Sp +Sp) represents the interaction of a slow wave with 
an Alfven wave travelling in one direction along the mag- 
netic field to produce an Alfven wave travelling in the op- 
posite direction. This generation of "back-scattered" Alfven 
waves does not produce waves at fe, = and thus does not 
contribute to AAA interactions or the associated perpendicu- 
lar cascade of Alfven-wave energy. Although the (hypothet- 
ical) conversion of energy into A+ energy would violate 
cross-helicity conservation in incompressible MHD, the gen- 
eration of back-scattered Alfven waves in compressible MHD 
does conserve cross helicity when one takes into account the 
associated change in the average flow velocity (v,). 

The d{k— q) term in equation ^ that does not contain 
represents the generation of slow waves by fast waves. If 
the wave fields are viewed as the sum of wave quanta, each 
of energy ft|(Bi:| and momentum hk, then this b{k — q) term 
represents the process f ^ f + s, i.e., a fast-wave decay- 
ing into a slow wave and a new fast-wave. This 8{k ~ q) 
term conserves the total number of fast-wave quanta Nf = 
j cPk [Fj^ / {h\(iy^ but decreases the fast- wave energy £/ = 

/ (fikFk, and thus causes an inverse cascade of fast-wave 
quanta to smaller frequency, i.e., a decrease in the average 
fast-wave frequency o5/ = Ef/hNf. The energy drained from 
fast waves is transferred to slow waves [through the 5{p — q) 
term in equation (|7]i], which are then rapidly damped. 



The term {ti/Ava) S crpcrqh{k 



q)^{qz 



kz)Ak\A^{d/dqz){q,A^) in equation ^ is denoted 
and the corresponding term in the equation for dA^ /dt 
is denoted . These terms represent the generation 
of slow waves by the interaction of oppositely di- 
rected Alfven waves, i.e. a* ^ + s. Upon defining 
E"^ ~ J dkxdkyAf and Qk = Tik^E^E'j^ /va, one can show that 
Hk = Jdk,dky{I+ +i;) = -Qk + {d/dk,){L,Qk), where the 
{d /dkz){kzQk) term represents a flux of Alfven-wave energy 
to smaller (inverse cascade). The energy drained from 
Alfven waves via the —Qk term in Hk is transferred to slow 
waves through the 5(/7- — q-) term in equation which then 
undergo rapid ion Landau damping. |18] This mechanism 
for transferring Alfven-wave energy to the ions is weak for 
"quasi-2D" fluctuations with \kz\ <C k± because of the factor 
of fc? in There are additional terms containing in the 
equation for 8A J /dt that result in the transfer of Alfven-wave 
energy to larger \kz\ [19], but these terms are higher order 
than e^^ and are neglected in this Letter since colHsionless 
damping keeps small at nonzero k^. 

Equation dHJ can be modified to allow for the possibility of 
strong Alfven-wave turbulence at small by replacing the 
AAA term [0= 5(?z)] in equation ([S) with the advection and 



diffusion terms on the right-hand side of Eq. (15) of mul- 
tiplied by a factor of 2 to convert to the normalization of Af 
used in this Letter Similar generalizations are possible for the 
"phase-mixing" and "passive-scalar mixing" terms. 

The interplay between the various nonlinear processes de- 
scribed in this Letter depends upon the value of p as well as 
the amplitudes and anisotropics of the different wave types at 
the forcing scale or "outer scale." For example, greater excita- 
tion of Alfven waves with > k± and fast-waves strength- 
ens the parallel cascade. (Alfven waves at k± ^ \kz\ cause 
only a weak secondary excitation of the fast waves and Alfven 
waves with > k± that participate in the parallel cascade.) 
On the other hand, the perpendicular cascade is strengthened 
by increasing the excitation at k, = of S^, entropy waves, 
and both and A^ . A stronger perpendicular cascade then 
weakens the parallel cascade by draining energy out of the 
"quasi-parallel" region of A;-space in which |^-| > k \ , redu cing 
the amount of wave energy that reaches very large \kz\. jlln 
Numerical solutions to equations O through (|9]l will be useful 
for describing turbulence in settings such as the solar corona, 
solar flares, and Earth's magnetosphere. 
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